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Abstract 

A fully self-consistent microscopic framework for evaluation of nuclear weak-interaction rates 
at finite temperature is introduced, based on Skyrme functionals. The single-nucleon basis and 
the corresponding thermal occupation factors of the initial nuclear state are determined in the 
finite-temperature Skyrme Hartree-Fock model, and charge-exchange transitions to excited states 
are computed using the finite-temperature RPA. Effective interactions are implemented self- 
consistently: both the finite-temperature single-nucleon Hartree-Fock equations and the matrix 
equations of RPA are based on the same Skyrme energy density functional. Using a representative 
set of Skyrme functionals, the model is applied in the calculation of stellar electron-capture cross 
sections for selected nuclei in the iron mass group and for neutron-rich Ge isotopes. 

PACS numbers: 21.60.Jz, 23.40.Bw, 23.40.Hc, 26.50.+X 
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I. INTRODUCTION 



Nuclear weak-interaction processes play a crucial role in the late stages of evolution of 
a massive star and in the pre-supernova stellar collapse The core of a massive 

star at the end of hydrostatic burning is stabilized by electron degeneracy pressure, as long 
as its mass does not exceed the Chandrasekhar limit of about 1.44 solar masses. If this 
mass limit is exceeded, electron pressure can no longer stabilize the core and it collapses. 
The dynamics of this process depends on the core entropy and lepton-to-baryon ratio ji], 
which are essentially determined by weak-interaction processes: the nuclear /5-decay 

(Z, iV) (Z + 1, iV - 1) + e- + z/, , (1) 

and electron capture: 

e- + {Z,N) ^Ue + {Z -1,N + 1) . (2) 

The latter reduces the number of electrons available for pressure support, whereas beta de- 
cay acts in the opposite direction. At low matter densities p < 10^^ g cm~'^ neutrinos escape 
from the star, carrying away energy, i.e. cooling the stellar core, and keeping its entropy 
low. For initial values of ~ 0.5, /^-decays are impeded by the presence of electrons which 
reduce the available phase space for decay, but become competitive when the composition 
of nuclei in the core gets more neutron-rich. In the early stage of the collapse, for densi- 
ties lower than a few 10^'^ g cm~'^, the electron chemical potential is of the same order of 
magnitude as the nuclear Q value, and the electron-capture cross sections are sensitive to 
the details of the Gamow- Teller GT'*' strength distributions in daughter nuclei. At these 
densities and temperatures between 300 keV and 800 keV, electrons are captured on nuclei 
with mass number A < 60. For higher densities and temperatures T ^ 1 MeV, electron 
capture occurs on heavier nuclei A > 65. Under these conditions, however, the electron 
chemical potential is significantly larger than the nuclear Q- value, and the capture rates are 
dominantly determined by the centroid and the total GT+ strength 

The first standard tabulation of nuclear weak-interaction rates for astrophysical applica- 
tions was that of Fuller, Fowler and Newman [b^ . It was based on the independent particle 
model, but used experimental information whenever available. The tables included rates for 
electron capture, positron capture, /3-decay, and positron emission for relevant nuclei in the 
mass range 21 < A < 60. Based on data which in the meantime became available on GT"*" 
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strength distributions, and using large-scale shell-model diagonalization in the complete pf- 
shell, these rates have been improved, and rates for electron and positron captures, and for 
/?"'" and (3~ decays have been cornputed for stellar conditions and for more than 100 nuclei 
in the mass range A = 45 — 65 ja, H]- Using the improved weak-interaction rates, models 
for pre-supernova evolution of massive stars were examined in Ref . Q| , and it was concluded 
that the resulting changes in the initial value of Yg and iron core mass could have important 
consequences for nucleosynthesis and the supernova explosion mechanism. Detailed calcula- 
tions of stellar weak-interaction rates in the iron mass region have also been carried out with 
the shell model Monte Carlo (SMMC) approach 91]. The advantage of this approach is that 
it treats nuclear temperature exactly, and can even include larger model spaces. There are 
limitations, however, in applying the SMMC to odd-A and odd-odd nuclei at low temper- 
atures. In addition, the SMMC yields only an averaged GT strength distribution, whereas 
the diagonalization shell-model approach allows for detailed spectroscopy. 

At higher densities and core temperatures T ^ 1 MeV, the excitation energy of a nucleus 
with mass number A ^ 80 is much larger than the energy gap between the pf- and sdg- 
shells. Weak-interaction rates for nuclei beyond the p/-shell cannot yet be systematically 
evaluated with large-scale diagonalization shell-model calculations, because of huge config- 
uration spaces and the lack of a reliable effective interaction in this mass region. Thus in 



Ref. 



101 ] a hybrid model has been introduced in which the nucleus is described as a Slater 



determinant with temperature-dependent occupation numbers, determined with shell-model 
Monte Carlo (SMMC) calculations. In the second step the electron capture rates are com- 
puted from GT"*" strength distributions calculated with the random-phase approximation 
(RPA) built on top of the temperature-dependent Slater determinant. The SMMC/RPA 
hybrid model was used to calculate electron capture rates on nuclei with mass numbers 
A = 65 — 112, at temperatures and densities characteristic for core collapse [l])]. It was 
shown that these rates are so large that electron capture on nuclei dominates over capture 
on free protons, and this leads to significant changes in the hydrodynamics of core collapse 
and bounce jlll, Q • 

The latest theoretical and computational advances in modeling the nuclear physics input 
for astrophysical applications have highlighted the need for fully microscopic global predic- 
tions for the nuclear ingredients. This is especially important when considering neutron-rich 
nuclei far from the line of /5-stability, for which data on ground-state properties and exci- 
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tations are not available. The basic advantages of the shell model is the ability to describe 
simultaneously all spectroscopic properties of low-lying states for a large domain of nuclei, 
and the use of effective interactions that can be related to two- and three-nucleon bare forces. 
On the other hand, since effective interactions strongly depend on the choice of active shells 
and truncation schemes, there is no universal shell-model interaction that can be used for all 
nuclei. Moreover, because single-particle energies and a large number of two-body matrix 
elements have to be adjusted to data, extrapolations to nuclei far from stability are not 
expected to be very reliable. Medium-heavy and heavy nuclei with very large valence spaces 
require calculations with matrix dimensions that are far beyond the limits of current shell 
model variants. Properties of heavy nuclei with a large number of active valence nucleons 
are therefore best described in the framework of nuclear energy density functionals (NEDF). 
At present NEDF's provide the most complete description of ground-state properties and 
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14( 1 . At the level of practical applica- 



coUective excitations over the whole nuclide chart 
tions the NEDF framework is realized in terms of self-consistent mean-field models. With 
a small set of universal parameters adjusted to data, this approach has achieved a high 
level of accuracy in the description of structure properties over the whole chart of nuclides, 
from relatively light systems to superheavy nuclei, and from the valley of /5-stability to the 
particle drip-lines. 

In this work we introduce a fully self-consistent microscopic framework for calculation 
of weak-interaction rates at finite temperature, based on Skyrme functionals. The single- 
nucleon basis and the corresponding thermal occupation factors of the initial nuclear state 
are determined in the finite-temperature Skyrme Hartree-Fock model, and charge-exchange 
transitions to excited states are computed using the finite-temperature RPA. Effective in- 
teractions are implemented self-consistently, i.e. both the finite-temperature single-nucleon 
Hartree-Fock equations and the matrix equations of RPA are based on the same Skyrme 
energy density functional. The advantage of this approach over shell-model calculations or 
hybrid models, is that a particular finite-temperature Hartree-Fock plus RPA model, i.e. 
completely determined by the choice of a Skyrme functional, can be extended over arbi- 
trary mass regions of the nuclide chart, without additional assumptions or adjustment of 
parameters, as for instance single-particle energies, to transitions within specific shells. In a 
simple RPA, of course, correlations are described only on the one-particle - one-hole level, 
and therefore one cannot expect the model to reproduce the details of the fragmentation 
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of GT strength distributions. This can only be accomphshed in the shell-model approach 
which includes higher-order correlations. In general, however, the RPA reproduces the cen- 
troid of strength distributions and the total GT strength. For electron capture, particularly, 
the RPA is an appropriate tool for the evaluation of cross sections for capture on nuclei in 
conditions where the electron chemical potentials are larger than the characteristic nuclear 
Q- values j^. 

Rather than evaluating and tabulating weak-interaction rates for hundreds of nuclei al- 
ready at this stage, in the present work we perform illustrative calculations of electron- 
capture cross sections for selected nuclei in the iron mass group and for neutron-rich Ge 
isotopes, and compare results with those obtained with the SMMC approach p] and the 
hybrid SMMC/RPA model 10|, respectively. Calculations are performed for a represen- 
tative set of Skyrme functionals, and this provides an estimate of the range of theoretical 
uncertainty inherent in the Skyrme Hartree-Fock plus RPA approach. 

The framework of finite-temperature Skyrme Hartree-Fock plus charge-exchange RPA, 
and the formalism for calculating cross sections for electron capture, are introduced in Sec- 
tion [Tll Electron capture on iron group nuclei {A ^ 45 — 65) is considered in Section UTTl and 
cross sections for electron capture on neutron-rich Ge isotopes are evaluated in Section IIV[ 
Section |V] summarizes the results of the present investigations and ends with an outlook for 
future studies. 



II. CALCULATION OF ELECTRON CAPTURE CROSS SECTIONS WITH 
FINITE-TEMPERATURE SKYRME HARTREE-FOCK PLUS RPA 

A. Charge-exchange RPA at finite temperature 

Throughout pre-supernova evolution, electron capture on nuclei proceeds at finite temper- 
ature. In addition to capture on p/-shell nuclei, this process also takes place on neutron-rich 
nuclei with protons in the p/-shell and neutron number > 40. Finite temperature effects 
and correlations unblock Gamow- Teller transitions that are forbidden at zero temperature. 

In the present analysis we employ the fully self-consistent finite temperature charge- 
exchange random phase approximation (FTRPA), formulated in the single-nucleon basis of 
the Skyrme Hartree-Fock model at finite temperature (FTSHF). Effective interactions are 
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implemented self-consistently, i.e. both the FTSHF equations and the matrix equations of 
FTRPA are based on the same Skyrme energy density functional. For a description of open- 
shell nuclei it is also necessary to include a consistent treatment of pairing correlations like, 



for instance, in the finite temperature HFB+QRPA framework 
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16j. However, in nuclei 

the phase transition from a superfluid to normal state occurs at temperatures T ^ 0.5 — 1 
MeV 17, [is], [3], whereas for temperatures above T ^ 4 MeV contributions from states in 
the continuum become large, and additional subtraction schemes have to be implemented to 
remove the contributions of the external nucleon gas 1201. In this work we consider a range 

n 

of temperatures relevant for the stellar electron capture process: T = 0.5 — 1.5 MeV [101], for 
which the FTSHF plus FTRPA should provide an accurate description of the Gamow- Teller 
and forbidden transitions. 

The finite temperature HF framework 



20 



2l[ | has been successfully used in nuclear struc- 



ture calculations for many years. In the case of a Skyrme functional the finite temperature 
HF equations have the same form as at zero temperature, but the density reads 



(3) 



where, in addition to bilinear products of single-nucleon HF wave functions 0q,, the con- 
tribution of each single-nucleon state is determined by the Fermi-Dirac (FD) occupation 
factors 

f ce — ~, ; TTTPrVTrTTTT • (^) 



1 + g(l/fcT)(e„-M) ■ 

ea are the single-nucleon energies, and the chemical potential /i is determined by the conser- 
vation of the number of nucleons X)« fa = ^- In this sense the implementation of the finite 
temperature formalism is very similar to the treatment of pairing correlations in the BCS 
framework, with Fermi-Dirac factors replacing the occupancies f^. The Fermi-Dirac factors 
determine contributions from individual single-nucleon states to other types of densities as 
well, e.g. the kinetic energy density and the spin-orbit density. 

A detailed derivation of the (Q)RPA formalism at finite temperature can be found in 



Refs. 
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22 



23 



2J], and the first self-consistent extension to open-shell nuclei (finite tem- 
perature QRPA) has been recently reported in Ref. 1^. Finite-temperature linear response 



theory and RPA have been success: 
decay of hot nuclei 
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201, 



22 



23 



uL 
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y used in numerous studies of giant resonances and 

23, y , y be 
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3l| . The FTRPA represents 



the small amplitude limit of the time-dependent mean-field theory at finite temperature. 



Starting from the response of a time-dependent density matrix p{t) to a harmonic external 
field f{t) [3^, the equation of motion for the density operator reads 



idtp=[h[p] + f{t),p]. 



In the small amplitude limit the density matrix is expanded to linear order 



(5) 



m = p'+ 5p{t) 



where 



5p{t) = 5p(+)e-^^* + 6p^^^e 



^(^) +iujt 



and denotes the stationary ground-state density 



(6) 



(7) 



(8) 



and includes the thermal occupation factors of single-particle states fk- For the analogous 
expansion of the Hamiltonian operator h(t) = h^^^ + 5h{t), the linearized equation of motion 
reads 



+ 



5h 



6p 



5p,p 



(0) 



(9) 



ih dtSp= [U^\6p 

Taking the matrix elements of this equation between the states and |0), we obtain 

(ca - ep) 5p^^^ + ^(/^ - fs)v^^-^s5py + {fs - f^)v^-s-p^5p^~^ , (10) 
for terms multiplying the factors e~*'^* and e"'"*'^*, respectively. 



-huj6pif3^ 



^h^ap = J2ifs - fl)'VaSp^Sp^^^ + (/^ - fs)v^^-ps5p\s^ (11) 
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and analogously for 6h[^p . These relations are consistent with the definition of the HF mean 
field at finite temperature, and express the fact that thermal occupancies determine the 
way density fluctuations affect the mean field. Equations fllUj) are consistent with those 
defined in Ref. |22|. If the first (second) equation of the set ffTUl) is multiplied by fp — fa 
{fa — f/3), with the definition F^/j = Sp^^J and = ^P^aJ ^ then equations (ITU!) take the 
. The set of equations (fTOj) is also consistent with the derivation of finite 



form as in Ref. 



temperature QRPA in Ref. 



15| . in the limit of vanishing pairing correlations. However, in 



15l | a different definition of tlie RPA amplitudes is introduced, wliicli involves quantities like 



y fa ~ fi3- This requires special care in the proton- neutron case because no simple condition 
can be imposed which would guarantee that the quantity under the square root is positive. 

The finite temperature forward- and backward-going amplitudes can be related to the 
corresponding zero-temperature amplitudes {X and Y) through the following relations 



The charge-exchange RPA matrices are composed of matrix elements of the residual in- 
teraction V, as well as certain combinations of thermal occupation factors fk- Because of 
finite temperature, the configuration space includes particle-hole (ph), particle-particle (pp), 
and hole-hole [hh) proton-neutron pairs. The residual interaction is derived from a Skyrme 
energy density functional, and single-particle occupation factors at finite temperature are 
included in a consistent way both in the FTSHF and FTRPA. The same interaction is used 
both in the FTSHF equations that determine the single-nucleon basis, and in the matrix 
equations of the FTRPA. The full set of FTRPA equations is solved by diagonalization. 
The result are excitation energies, and the corresponding forward- and backward-going am- 
plitudes that are used to evaluate the transition strength for a given multipole operator. 
As an illustrative example, in Figs. [1] and [2] we display the temperature dependence of the 
occupations of the proton and neutron orbitals fj/2, P3/2, h/2, P1/2 and gQ/2 in '^''Ge, calcu- 
lated in the finite-temperature Skyrme HF model. The self-consistent calculation with the 
SGII effective interaction [sS^ corresponds to zero temperature, T = 0.5 MeV and T = 1.3 
MeV. At zero temperature the proton orbitals f7/2 and P3/2 are fully occupied, whereas f5/2, 
Pi/2 and g9/2 are empty. By increasing the temperature protons are mostly promoted from 
the P3/2 into the f5/2 orbital and, to a lesser extent, into the pi/2 orbital. Correspondingly 
the occupation percentage of P3/2 is reduced, whereas those of the f5/2 and pi/2 orbitals are 
enhanced. In the temperature intervals considered here, the occupations of the f7/2 and g9/2 
orbitals do not change significantly. At zero temperature neutrons fully occupy the f7/2, P3/2, 
f5/2, and pi/2 orbitals, and there are two neutrons in the g9/2 orbital. The main effect of 
increasing temperature is to promote neutrons from P3/2, f5/2, and pi/2 into g9/2. However, 
note that for protons already at T = 0.5 MeV the calculation predicts a pronounced effect 




— ^a/sffsi^ " fa) + YiSafai^ — fp) , 



(12) 



and 



¥a~J = Yapfp{l - fa) + XpJ^(l - fp) . 



(13) 
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FIG. 1: (Color online) Occupation percentage of the proton orbitals £7/21 P3/2) ^5/2? P1/2 ^^'^ S9/2 
in ^^Ge, calculated in the finite-temperature Skyrme HF model with the SGII interaction, at zero 
temperature, T = 0.5 MeV and T = 1.3 MeV. 
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FIG. 2: (Color online) Occupation percentage of the neutron orbitals f7/2, P3/2) h/2i P1/2 gg/2 
in ^^Ge, calculated in the finite-temperature Skyrme HF model with the SGII interaction, at zero 
temperature, T = 0.5 MeV and T = 1.3 MeV. 
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FIG. 3: (Color online) Temperature dependence of the Gamow- Teller GT~ and GT"*" strength 
distributions in ^^Ge, calculated with the finite-temperature proton-neutron RPA model based on 
the Skyrme SGII interaction. 

on the occupation of orbitals close to the Fermi surface, whereas the occupation of neutron 
orbitals is significantly modified only at the higher temperature of T = 1.3 MeV. The re- 
sults shown in Figs. [T] and [2] can be compared with thermal occupation numbers obtained 
from canonical shell model Monte Carlo (SMMC) calculated in Ref. lOj, with Woods-Saxon 
single-particle energies and a pairing plus quadrupole residual interaction. The temperature 
dependence of occupation numbers predicted by the two models is similar, with the SMMC 
results showing a more pronounced effect on occupation numbers already at T = 0.5 MeV, 
especially for neutron orbitals. This can be attributed to a smaller energy gap between 
the pi/2 and g9/2 single-neutron orbitals used in the SMMC calculation, and to additional 
correlations in the ground state that are not taken into account in the simple Skyrme HF 
model. 

In Fig. [3] we display the temperature dependence of the corresponding Gamow- Teller 
GT^ and GT+ strength distributions in ^^Ge, calculated with the finite-temperature proton- 
neutron RPA model based on the Skyrme SGII interaction. In the GT^ direction a neutron 
is changed into a proton, as in /3-decay. At low temperature the GT~ mode corresponds 
to a coherent superposition of J'^ = charge-exchange proton-particle - neutron-hole 
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transitions. The GT operator reads 



t$t=j:^^±- (14) 



In addition to the high-energy GT^ resonance at ~ 14 MeV, a collective superposition 
of direct spin-flip (j = I + ^ j = I — \) transitions, the response function displays a 
concentration of strength in the low-energy tail. The transitions in the low-energy region 
correspond to core-polarization (j = / ± | ^ j = / ± and back spin-flip {j = I — ^ 
j = / + |) neutron- hole - proton-particle excitations. The GT^ transitions are allowed 
at T = 0, and the distribution displays only a weak temperature dependence. On the 
other hand, for ^"^Ge the GT"*" transitions, in which a proton is changed into a neutron, 
are forbidden at zero temperature because the relevant neutron orbitals are fully occupied. 
Note that additional ground-state correlations, e.g. pairing, are not taken into account in 
the Skyrme Hartree-Fock plus RPA model. Even at T = 0.5 MeV the occupation factors 
of neutron orbitals remain unchanged (cf. Fig. [2]) and the calculation does not predict any 
low-energy GT+ transition. Only at higher temperatures above T = 1.0 MeV the low-energy 
GT+ transitions are thermally unblocked. We have verified that at each temperature the 



Gamow- Teller strength distribution satisfies the Ikeda sum rule 34 1 



{Sgt - S^t) = m - Z) , (15) 



where Sqj, denotes the total sum of Gamow- Teller strength for the GT^ transitions. 



B. Cross section for electron capture 



The electron capture on a nucleus [Z, N) 



e- + (Z, N) -^Ue + iZ -1,N+ ly , 



(16) 



presents a simple semi-leptonic reaction that proceeds via the charged current of the weak 
interaction. The theoretical analysis of these processes necessitates the description of the 
weak interaction between leptons and nucleons, as well as the wave functions of the initial 
and final nuclear states. Detailed expressions for the reaction rates and the transition matrix 



elements can be found in Refs. 



35 



36 



37l |. The electron-nucleus reaction cross section for 
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a transition between the states \i) and |/) reads 

l^-H, (17) 

^ ' lepton spins ' Mi Mf 

where V denotes the quantization volume, and Ei, is the energy of the outgoing electron 
neutrino. The Hamiltonian Hw of the weak interaction is expressed in the standard current- 
current form, i.e. in terms of the nucleon J\{x) and lepton j\{x) currents 

Hw = / d'x Jx{x)j\x) , (18) 

where G is the weak coupling constant, and the resulting transition matrix element reads 

mw\^ = -^h J ^ e-^''-{f\j\x)\z) , (19) 

where the four-momentum transfer is g = (go, q), and the multipole expansion of the leptonic 
matrix element lxe~'^'^'^ determines the operator structure for the nuclear transition matrix 



elements 



35 



36 



33] . The expression for the electron capture cross sections is given by 



da _ G^cos^e^ F{Z, E^) 
dn ~ 2^ (2 Ji + 1) 

X 



I J>1 

-2q-(^-/3)i?e(J;||r/^^||J.)(J;||r/^||J,,)*} 



+ EW(£^e,^.){(l-'>-/3 + 2(i>-q)(/3-q)(J/||i:j||Ji)|2 + (l + i>-/3)(J/||^j||Ji)|2 

J>0 

-2q{i^ + P)Re{Jf\\jCjm{Jf\\MjUy}'^ , (20) 

where the momentum transfer q = u ~ k is defined as the difference between neutrino and 
electron momenta, q and u are the corresponding unit vectors, and f3 = k/E^. The energies 
of the incoming electron and outgoing neutrino are denoted by Eg and E^, respectively. The 
Fermi function F{Z, E^) corrects the cross section for the distortion of the electron wave 
function by the Coulomb field of the nucleus 



'^(^") - (Trim • 

with the phase-space factor (1 + E^/Mx)~^ accounting for the nuclear recoil, and Mj- is the 
mass of the target nucleus. The nuclear transition matrix elements between the initial state 
I Jj) and final state \Jf), correspond to the charge Aij, longitudinal Cj, transverse electric 
Tj^^, and transverse magnetic Tf^^^ multipole operators: 
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the Coulomb operator 



Mjm{x) = F^Mf\x) 
the longitudinal operator 



K 



rriN 



F^nfix) + -{Fa - m,Fp)T.f\x) 



(22) 



K 



(23) 



the transverse electric operator 



K 



rriN 



+ iFjj:'f{x) 



(24) 



and the transverse magnetic operator 



F,-Af(a.)--^^E'f(a.) 



+ FA^'f{x 



(25) 



where all the form factors are functions of g^, and k = \q\. The operators M, fi, A, A', 
S, S', and S" are expressed in terms of spherical Bessel functions, spherical harmonics, and 
vector spherical harmonic s [35i] . By assuming conserved vector current (CVC), the standard 
set of form factors reads |39|: 



W)= 1 + 

/i^(g2) = 4.706 
FA{q^) = -1.262 
Fpiq') 



840 MeV 



1 + 



1 + 



840 MeV J 

q 

1032 MeV 



-2 



2, 2mNFA{q^^ 



q2 _^ ^2 



(26) 

(27) 
(28) 
(29) 



The cross sections for electron capture are evaluated from Eq. fl20|) . with transition matrix 
elements between the initial and final states determined in a self-consistent microscopic 
framework based on the (finite temperature) Skyrme HF model for the nuclear ground 
state, and excited states are calculated using the corresponding (finite temperature) RPA. 
For each transition operator Oj the matrix elements between the initial state of the even- 
even {Z, N) target nucleus and the final state in the corresponding (Z — 1, 1) nucleus are 

13 



expressed in terms of single-particle matrix elements between the single-nucleon states, and 

'^p^ and 5p^~^' 



the corresponding (finite temperature) RPA amplitudes 5p~^)'^ and Sp^^J"^ (cf. Sec. Ill Ap 



(^/l|o,||J.) = Ei^Wom (5pS - ^pS ) ■ (30) 

The energy of the outgoing neutrino is determined by the conservation relation: 

E, = Ee-Q + Ei-Ef , (31) 

which includes the difference between the final and initial nuclear states. The Q-value plays 
a particularly important role in the calculation of electron capture rates. Namely, the energy 
that is available to excite states in the daughter nucleus depends on whether electron capture 
on a specific target nucleus releases energy (Q < 0), or requires an additional external input 
Q > 0). In the present calculation the Q-value is determined from the experimental masses 



40(1 : Q = Mf—Mi, where Mjj are the masses of the parent and daughter nuclei, respectively. 

The nuclei that will be considered in this work contribute to stellar electron capture rates 
in the temperature interval T ^ 0.5 — 1.5 MeV [10]. The expression for the total cross 
section for electron capture on a nucleus {Z, N) at temperature T reads 

.iE., D ^ g E FiZ, E.)iB±^£^j: iE.-Q^E.- E,fl^^ , (32) 

where Oj is the generic notation for the charge Aij, longitudinal Cj, transverse electric Tf'^, 
and transverse magnetic Tf^^^ multipole operators. The sum over initial states includes a 
thermal average of levels, with the corresponding partition function G{Z, A,T). The finite 
temperature induces the thermal population of excited states in the parent nucleus. Each 
of these states \i) is connected by the multipole operators to many levels |/) in the daugh- 
ter nucleus. The calculation of all possible transitions is computationally prohibitive, and 



there 





ore the evaluation of the total cross section for electron capture is usually simplified 



10 



41| by adopting the Brink hypothesis, i.e. by assuming that the strength distribu- 
tion of the multipole operators in the daughter nucleus is the same for all initial states and 
shifted by the excitation energy of the initial state. By using this approximation, the sum 
over final states becomes independent of the initial state and the sum over the Boltzmann 
weights cancels the partition function. The Brink hypothesis is a valid approximation when 
the temperature and density are high enough so that many states contribute, and variations 



14 



in the low-energy transition strength cancel out. As it was done in the calculation of stellar 



electron capture on neutron-rich germanium isotopes in Ref. [10|, we apply the Brink hy- 
pothesis to the initial state which represents the thermal average of many-body states in the 
parent nucleus at temperature T. This thermally averaged initial state is approximated by 
the finite temperature Skyrme Hartree-Fock ground-state Slater determinant with Fermi- 
Dirac thermal occupation factors. With this approximation the final expression for the total 
electron capture cross section at temperature T reads 

cr(Ee, T) = ^FiZ, E,) J2{E,-Q- ujff ^ Sj{ujf, T) , (33) 

/ J 

where 0;^ is the excitation energy in the daughter nucleus, and Sj is the discrete finite- 
temperature RPA response for the multipole operator Oj. 

III. ELECTRON CAPTURE ON IRON GROUP NUCLEI 

As our first illustrative example and application of the model, we consider electron capture 
on iron group nuclei {A ^ 45 — 65). The calculated electron capture cross sections and rates 
for nuclei in this mass range are essential for modeling the initial phase of stellar core-collapse 
and supernova explosion nay , |4l|]. In the presupernova collapse electron capture on pf- 



shell nuclei proceeds at temperatures between 300 keV and 800 keV. Detailed calculations of 
stellar weak-interaction rates in the iron mass region have been carried out in the framework 
of the interacting shell model j^. Both the shell model Monte Carlo (SMMC) ^, and large- 



scale shell-model diagonalization [6|, [?], were used to calculate electron capture and /5-decay 
rates in the A ^ 45 — 65 mass region. The SMMC results [3] were in fact superseded in 
Ref. [g], by the weak-interaction rates obtained using large-scale shell-model diagonalization 
in the complete p/-shell. 

In this section we calculate electron capture cross sections for selected nuclei in the iron 
mass region, and compare the results with those of Ref. j^, where the SMMC was used to 
calculate Gamow- Teller GT'*' strength distributions (in this direction a proton is changed 
into a neutron), and these distributions were then used to compute the electron-capture 
cross sections and rates in the zero-momentum transfer limit, as a function of the incident 
energy of the electron. The SMMC calculations solve the full shell-model problem for the 
GT+ strength distributions in the Ohu /p-shell space using a realistic residual interaction. 
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FIG. 4: (Color online) Electron capture cross sections for the reaction ^^Fe(e~, fe)^^Mn at T = 0.5 
MeV, as functions of the incident electron energy E^. The Skyrme HF+RPA results are compared 
with cross sections calculated from the SMMC Gamow- Teller strength distributions ^. 



a 



In the calculation of Ref. [9| the KB3 residual interaction 42] was used, and the quenching 
of the total GT strength was taken into account by renormalizing the GT transition matrix 
elements by the constant factor 0.8. In the present analysis the cross sections for electron 
capture are evaluated from Eq. ( |33l) . with transition matrix elements between the initial 
and final states determined in the self-consistent microscopic framework based on the finite 
temperature Skyrme HF model for the nuclear ground state, and excited states are calculated 
using the corresponding finite temperature charge-exchange RPA. 

In Fig. m we display the calculated cross sections for the reaction ^^Fe(e~, t'e)^^Mn as 
functions of the incident electron energy Eg. Calculations are performed at temperature 
T=0.5 MeV, and the Skyrme HF+RPA results are compared with cross sections calculated 
from SMMC Gamow- Teller strength distributions j^. As the calculation of cross sections 
in Ref. 9j corresponds to the zero-momentum transfer limit and only includes the Gamow- 
Teller operator, for the sake of comparison we have also limited the sum in Eq. fl20|) to 
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the l"*" channel, i.e. only transitions to l"*" excited states are taken into account. At low 
momentum transfer, the allowed Gamow- Teller transitions dominate the electron capture 
process on p/-shell nuclei. Note, however, that in the calculation of Dean et al. only 
the Ohuj Gamow- Teller transition strength is considered, rather than the total strength in 
the l"*" channel. The reduction of the axial-vector coupling constant from its free-nucleon 
value gA = 1.262 to = 1-0 (cf. Eq. (!28|) ). is equivalent to the renormalization of the 



GT matrix elements in Ref. 
weak-interaction rates in t 
diagonalization approach 



9| by the constant factor 0.8. In the shell-model studies of 
le A ~ 45 — 65 mass region, both the SMMC and shell-model 
^, have used the KB3 residual interaction, which is well suited 
for full Ohuj calculations throughout the lower p/-shell region. However, to calculate the 
weak-interaction rates in the entire mass range A ^ 45 — 65, the original KB3 interaction 
had to be modified by including a number of monopole corrections in order to reproduce the 
GT strength distributions and half-lives. 

In the present analysis the cross sections, as functions of the incident electron energy, are 
compjuted for a representative set of Skyrme functionals: SGII [33], SkM* [43], SLy4|44|. 
SLy5[45|], and SkO' 461]. Over the last thirty years more than hundred different Skyrme 
parameterizations have been adjusted and analyzed, and it is often difficult to compare 
results obtained with different models, also because they include different subsets of terms 
from the most general functional. Since in this work we apply the microscopic approach 
based on Skyrme HF+RPA, calculations are performed using various Skyrme functionals. 
In principle this will provide an estimate of the range of theoretical uncertainty inherent 
in the present approach. The electron-capture cross sections in Fig. H] exhibit a sharp 
increase of several order of magnitude within the first few MeV above threshold, and this 
refiects the GT"*" distributions. For electron energy Ee > 10 MeV the calculated cross 
sections display a more gradual increase. A very similar energy dependence of the cross 
sections is calculated for the neighboring even-even parent nuclei ^^Ti and ^°Cr, in Figs. [5] 
and El respectively. At low energies all Skyrme HF+RPA cross sections are below the 
values based on SMMC calculations. This is especially pronounced in ^^Fe(e~, z/e)^^Mn, and 
much less in ^°Cr(e", z/e)^°V. Cross sections calculated at very low electron energies will be 
very sensitive to the discrete level structure of the Gamow- Teller transitions. These cross 
sections, however, are several orders of magnitude smaller than those for E^. > 10 MeV and, 
when folded with the electron fiux to calculate capture rates, the differences between values 
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FIG. 5: (Color online) Same as in Fig. [U but for the reaction ^^Ti(e , fe)^^Sc . 

predicted by various models in the low-energy interval will not have a pronounced effect 
on the electron capture rates j^. Note, however, that in general this will strongly depend 
on the matter density and temperature of the environment. More important could be the 
differences at higher electron energies > 10 MeV, for which the Skyrme HF+RPA model 
systematically predicts cross sections above the values based on the SMMC The reason for 
this systematic effect is most probably that SMMC calculations are carried out only in the 
Ohuj /p-shell space. Note also that the spread in the calculated cross sections at low energies 
is greatly reduced for higher incident energies and, above ^ 15 MeV, all Skyrme effective 
interactions used in the present HF+RPA calculation effectively predict a universal behavior 
of the total electron capture cross sections. 

In general, for heavier nuclei and higher electron incident energies, not only the l"*", but 
other multipole transitions will also contribute to the total cross section for electron capture 
(cf. the sums in Eq. ( |20l) ). In calculations for the iron group nuclei based on the shell- model, 
;he cross sections were calculated only from the Ohu Gamow- Teller strength distributions 
6|,l9|. In Fig. [7] we plot the electron capture cross section for the ^^Fe target nucleus at T=0.5 
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FIG. 6: (Color online) Same as in Fig. HJ but for the reaction ^'^Cr(e , i^e)^'^V. 

MeV, calculated with the FTHF+FTRPA using the SLy5 Skyrme effective interaction. In 
addition to the total cross section which includes multipole transitions J'^ = 0"^, 1^, and 
2^, contributions from the individual channels are shown in the plot, as functions of the 
incident electron energy Ee- In this case all the way up to ^ 30 MeV the total cross 
section is completely dominated by the l"*" channel, with contributions from other channels 
being orders of magnitude smaller. Only at very high electron energies contributions from 
other multipole transitions become sizeable. 

Finally, in Fig. [H] we illustrate the isotopic dependence of electron capture cross sections. 
Cross sections for even-even Ni target nuclei, i.e. for the reactions '^Ni(e~, h'e)^Co, {A = 56 — 
62) at T=0.5 MeV, are calculated in the FTHF+FTRPA with the Skyrme SLy5 interaction. 
At any given incident electron energy the calculated cross sections decrease systematically 
along the Ni isotope chain, because more and more neutron orbitals become occupied, and 
therefore not accessible for electron capture reactions. Since with the increase of the number 
of neutrons the Q value for electron capture increases, more energetic electrons are required 
for the capture reaction on neutron rich isotopes and, at any given electron energy, the total 
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FIG. 7: (Color online) Electron capture cross section for the ^^Fe target nucleus at T = 0.5 MeV, 
calculated with the FTHF+FTRPA using the SLy5 Skyrme effective interaction. In addition to 
the total cross section which includes multipole transitions J'^ = 0^, 1^, and 2^, contributions 
from the individual channels are shown in the plot as functions of the incident electron energy E^. 

cross sections are smaller. We note that the predicted isotopic dependence of electron capture 
cross sections in Ni nuclei is in qualitative agreement with the results of the SMMC-based 



study of Ref. 



9| 



IV. STELLAR ELECTRON CAPTURE ON NEUTRON-RICH GE ISOTOPES 

At higher densities and temperatures T 1 MeV during the collapse phase, electrons are 
also captured on heavier and more neutron-rich nuclei with protons in the p/-shell {Z < 40) 
and neutrons > 40. In a naive independent particle picture the Gamow- Teller transitions 
which, as it was shown in Sec. IIIH dominate electron capture in the p/-shell, are forbidden 
for nuclei with Z < AO and A^ > 40. However, as it has been demonstrated in several 
studies, GT transitions in these nuclei are unblocked by finite temperature excitations. A 
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FIG. 8: (Color online) Electron capture cross sections for the even-even Ni target nuclei (A=56-62) 
at r = 0.5 MeV, calculated with the FTHF+FTRPA using the SLy5 Skyrme effective interaction. 



very detailed study based on the random phase approximation 471], has shown that electron 
capture on nuclei with protons in the p/-shell and > 40 can compete with capture on free 
protons if forbidden transitions are taken into account, in addition to allowed ones. At high 
temperatures T ~ 1.5 MeV Gamow- Teller transitions are thermally unblocked as a result 
of the excitation of neutrons from the p/-shell into the (79/2 orbital. In Ref. [l^ a hybrid 
model has been introduced to calculate electron capture rates on neutron-rich nuclei in this 
mass region. In the hybrid model the temperature and configuration- mixing effects are 
taken into account with the Shell Model Monte Carlo (SMMC) method, and are described 
by partial occupation numbers for the various single-particle orbits. Using mean-field wave 
function with finite temperature occupation numbers determined from SMMC, the electron 
capture cross sections are calculated with an RPA approach. Both allowed GT and forbidden 
transitions are included in the calculation. The SMMC/RPA hybrid approach was applied 
to the even germanium isotopes ^^^''^Ge at typical core collapse temperatures T ~ 0.5 — 1.3 
MeV, and it was demonstrated that configuration mixing is strong enough to unblock the 
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FIG. 9: (Color online) Electron capture cross sections for the reaction ^^Ge(e^, i/e)^^Ga at T=0.5 
MeV, as functions of the incident electron energy Eg. The Skyrme HF+RPA results are compared 
with cross sections calculated with the hybrid SMMC/RPA model {lo| . 



Gamow- Teller transitions at all temperatures relevant to core-collapse supernovae pi| . The 
SMMC/RPA model was also used to calculate rates for electron capture on nuclei with mass 
numbers A = 65 — 112, at temperature and densities characteristic for core collapse jll |. 
It was shown that electron capture on nuclei dominates over capture on free protons, and 
simulations of core collapse have demonstrated that these capture rates produce a strong 
effect on the core collapse trajectory and the properties of the core at bounce. 

In the present analysis we apply the self-consistent finite-temperature Skyrme HF+RPA 
model in the calculation of stellar electron-capture cross sections on neutron-rich Ge nuclei. 



and compare the resu 
SMMC/RPA model 



ts with those obtained by Langanke, Kolbe and Dean, using the hybrid 



lol | . In Fig. [9] we plot the various Skyrme HF-I-RPA electron capture 
cross sections for the reaction ^^Ge(e~, z/e)''^Ga at T=0.5 MeV, as functions of the incident 
electron energy E^,, in comparison with the SMMC/RPA results. As in the calculation for 
the iron group nuclei in Sec. Illlt the set of Skyrme effective interactions includes: SGII 
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33|, SkM* |43|, SLy4|4^, SLy5|45|, and SkO' [46|. A given interaction is consistently used 



both in the finite temperature SHF equations that determine the single-nucleon basis, and 
in the matrix equations of the finite temperature RPA. The SMMC calculation of Ref. [l^ 
included the complete (p/5'9/2) shell-model space and used a pairing+quadrupole residual 
interaction with parameters adjusted for this mass region. The single-particle energies were 
adopted from the KB3 interaction, but the /5/2 orbital was artificially reduced by 1 MeV 
to simulate the effects of the ar component that is missing in the residual interaction. An 
energy splitting of 3 MeV between the (79/2 and the /5/2 orbitals was assumed. For the RPA 
calculation based on SMMC, the single-particle energies were taken from a Woods-Saxon 
parameterization, and the residual interaction is a finite-range G-matrix derived from the 



Bonn nucleon-nucleon potential 48|] . In the calculation of total cross sections, both models 
include the multipole transitions: J'^ = 0^, 1^, and 2^. In the SMMC/RPA calculation the 
GT strength is quenched by multiplying the GT transition matrix elements by the constant 
factor 0.7. In the present analysis the standard reduction of the axial vector coupling 
constant is employed = 1.262 ^ qa = 1-0, which corresponds to the quenching factor 
0.8. 

Similar to the case of iron group nuclei, the cross sections calculated with different Skyrme 
parameterizations display a spread of values of less than an order of magnitude at lower elec- 
tron energies. At higher incident energies the differences between values calculated with dif- 
ferent effective interactions are much smaller. In general the Skyrme FTHF-I-FTRPA results 
are in good agreement with the cross sections calculated in the hybrid SMMC/RPA model, 
especially at higher electron energies Ee > 20 MeV. At lower energies the SMMC/RPA cross 
sections are considerably above the results obtained in the present calculation. The Skyrme 
FTHF+FTRPA calculations have also been carried out at higher temperature: T = 1.3 
MeV (Fig. [TOj) . It is interesting to note that in this particular case the FTHF+FTRPA 
results are indeed very close to those obtained with the SMMC/RPA model. Note that 
in Ref. [l^ values of the calculated cross sections are only reported for incident electron 
energies i^e > 12 MeV. 

In Figs. [TT] and [T2] the electron capture cross sections are shown for the ^'^Ge and ''^Ge 
target nuclei at T = 1.3 MeV, respectively. The Skyrme FTHF+FTRPA results are in 
qualitative agreement with the values calculated in the SMMC/RPA model. For ^^Ge, in 
particular, the Skyrme interactions divide into two branches: SLy4, SLy5 and SkM* predict 
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FIG. 10: (Color online) Same as Fig. [U but for the temperature T = 1.3 MeV. 




FIG. 11: (Color online) Same as Fig. El but for ^^Ge at T = 1.3 MeV. 
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FIG. 12: (Color online) Same as Fig. El but for ^^Ge at T = 1.3 MeV. 

cross sections that are systematically larger than those obtained from the SMMC/RPA 
model, whereas cross sections calculated with SGII and SkO' are below the SMMC/RPA 
results for electron energies < 30 MeV. The Skyrme FTHF+FTRPA cross sections 
are systematically smaller than the values predicted by the hybrid SMMC/RPA model for 
the target nucleus ^^Ge. The isotopic dependence of the electron capture cross sections, 
illustrated in Fig. [8] for the even-even Ni nuclei, is also observed in Figs. [10] - [I2] for the 
Ge isotopes. With the increase of the neutron number the threshold for electron capture 
is shifted toward higher electron energies, reflecting the change in the Q-value. For a given 
electron incident energy, the total cross section is reduced with the increase of the number 
of neutrons. 

Finally, in Fig. [13] we illustrate the temperature dependence of electron capture on ''^Ge. 
The cross sections are calculated with the FTHF+FTRPA model at three temperatures: 
T = 0.5, 1.3, and 2.0 MeV, using the SLy5 parameterization. The notable increase in the 
calculated cross sections occurs between T = 0.5 and 1.3 MeV, and this corresponds to a 
significant thermal unblocking of the neutron P3/2, f5/2, P1/2 orbitals (cf. Figs. [1] and [2] for 
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FIG. 13: Electron capture cross sections for the Ge target nucleus at temperatures T = 0.5, 1.3, 
and 2.0 MeV, calculated with the FTHF+FTRPA using the SLy5 Skyrme effective interaction. 

the case of ^"^Ge). Because these orbitals are aheady unblocked at T = 1.3 MeV, a further 
increase in temperature to T = 2.0 MeV results only in a moderate enhancement of the 
electron-capture cross sections. 



V. SUMMARY AND CONCLUSIONS 



Recent advances in modeling nuclear structure phenomena have also had a strong impact 
on astrophysical applications. More and more often calculations of stellar nucleosynthesis, 
nuclear aspects of supernova collapse and explosion, and neutrino-induced reactions, are 
based on microscopic global predictions for the nuclear ingredients, rather than on simplified 
semi-empirical approaches. In general, the required nuclear input includes properties of 
hundreds of nuclei at and far from the line of /3-stability, including the characteristics of 
strong, electromagnetic and weak interaction processes. Many of these nuclei, especially on 
the neutron-rich side, are not accessible in experiments and, therefore, nuclear astrophysics 
calculations crucially depend on accurate theoretical predictions for the nuclear masses, 
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bulk properties, nuclear excitations, (n, 7) and (7,n) rates, a- and /3-decay half-lives, fission 
probabilities, electron and neutrino capture rates, etc. 

Improved microscopic stellar weak-interaction rates, evaluated with large-scale shell- 
model diagonalization and/or hybrid RPA models, have been employed in recent studies of 
pre-supernova evolution of massive stars, and it has been shown that the resulting changes in 
the lepton-to-baryon ratio and iron core mass lead to significant changes in the hydrodynam- 
ics of core collapse and the supernova explosion mechanism. These results have emphasized 
the need for accurate microscopic evaluations of nuclear weak-interaction rates, at densi- 
ties and temperatures characteristic for core collapse, that can be extended over arbitrary 
mass regions of the nuclide chart. In this work for the first time a self-consistent micro- 
scopic framework for calculation of weak-interaction rates at finite temperature has been 
introduced, based on Skyrme functionals. Single-nucleon levels, wave functions, and ther- 
mal occupation factors for the initial nuclear state are determined in the finite-temperature 
Skyrme Hartree-Fock model, and transitions to excited states are computed using the corre- 
sponding finite-temperature charge-exchange RPA. Effective interactions are implemented 
self-consistently, i.e. both the finite-temperature single-nucleon Hartree-Fock equations and 
the matrix equations of RPA are based on the same Skyrme energy density functional. 

The model has been employed in illustrative calculations for stellar electron capture on 
selected nuclei in the iron group mass region, and for neutron-rich isotopes of germanium. 
Electron-capture cross sections have been calculated as functions of the energy of the in- 
cident electron, for a representative set of Skyrme functionals. By using different Skyrme 
functionals one is able to estimate the range of theoretical uncertainty of the Hartree-Fock 
plus RPA approach. For the iron group nuclei, the results have been compared with those of 
Ref . ^ , where the SMMC was used to calculate Gamow- Teller GT"^ strength distributions, 
and electron-capture cross sections and rates were computed in the zero-momentum trans- 
fer limit. At low incident electron energies, at which the cross sections are sensitive to the 
discrete level structure of the Gamow- Teller transitions, all Skyrme HF+RPA cross sections 
are smaller than the values based on SMMC calculations. These cross sections, however, 
are very small and the differences between various models will not have a pronounced ef- 
fect on the calculated electron capture rates. More important could be the differences at 
higher electron energies i?e > 10 MeV, for which the Skyrme HF+RPA model systemati- 
cally predicts cross sections larger than the values evaluated with the SMMC. It has to be 
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emphasized that the RPA approach takes into account large configuration spaces so that for 
any multipole operator the whole sum rule is exhausted, whereas generally this is not the 
case in shell-model calculations. 

For electron capture on neutron-rich Ge nuclei, the finite-temperature Skyrme Hartree- 
Fock plus RPA cross sections have been analyzed in comparison with results obtained using 
the hybrid SMMC/RPA model [l^, in which the nucleus is described as a Slater deter- 
minant with thermal occupation numbers determined with the SMMC, and capture rates 
are computed using a charge-exchange RPA built on top of the temperature-dependent 
Slater determinant. In general, a very good agreement has been found between the Skyrme 
FTHF-I-FTRPA results and the cross sections calculated in the hybrid SMMC/RPA model, 
especially at higher electron energies E^. > 20 MeV, and higher temperatures T > 1 MeV. 
In all cases the two models predict a similar dependence of the cross sections on electron 
energy in the interval 12 < E^. < 30 MeV. There are, however, differences in the absolute 
values, especially at relatively low temperature (T = 0.5 MeV), electron energies E^ < 20 
MeV, and for heavier isotopes, e.g. ''^Ge. 

The results of the present study show that the finite-temperature Skyrme Hartree-Fock 
plus charge-exchange RPA framework provides a valuable universal tool for the evaluation 
of stellar weak-interaction rates. Based on a universal Skyrme energy density functional, in 
the sense that the same functional is used for all nuclei, this framework can be employed 
in studies of weak-interaction processes in different mass regions. At relevant incident elec- 
tron energies the absolute spread in the electron-capture cross sections, computed with 
a variety of Skyrme functionals, is less than an order of magnitude. The next step will 
be a more extensive calculation and tabulation of electron capture rates for nuclei in the 
mass range A ^ 60 — 80, which can be compared with modern semi-empirical estimates of 
weak-interaction rates for intermediate mass nuclei, calculated using available experimental 
information and simple estimates for the strength distributions and transition matrix ele- 
ments 491]. For open-shell nuclei at very low temperatures, and especially for calculation 
of /3-decay rates, the framework could be extended to include pairing correlations. In the 
zero-temperature limit, isobaric analog states and GT resonances in open-shell nuclei have 



recently been studied with the new^ 
RPA based on Skyrme functionals 



developed self-consistent quasiparticle charge-exchange 
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5l| . Another interesting extension of the model would 



be the use of a different class of nuclear energy density functionals, for instance relativistic 
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EDFs, in which case excited states could be calculated using the charge-exchange relativistic 
QRPA [521. 
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